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1.0  INTRODUCTION 


This  document  constitutes  the  Final  Report  on  Contract  N00014-90-C-0214/ 
P00001,  "Relative  Merits  of  Innovative  Imaging  Sensor  Systems,”  and  covers  the 
period  18  January  1991  through  12  April  1992.  The  purpose  of  the  work  was  to  con¬ 
duct  research  and  assist  efforts  in  ground-based  optical  identification  and  tracking  at 
the  Innovative  Science  and  Technology  Experimental  Facility  (ISTEF)  in  Florida.  The 
technology  areas  covered  included  signal  processing,  the  collection  of  target  signature 
data,  and  proof-of-concept  field  demonstrations  of  sensor  systems.  A  major  portion  of 
the  work  centered  on  the  design  of  a  laser  Doppler  field  experiment  which  would  be 
suitable  for  deployment  at  ISTEF,  a  site  which  is  uniquely  well  suited  for  obtaining  sig¬ 
nature  data  relevant  to  the  rocket  lift-off  episode.  Extensive  documentation  of  interim 
results  is  contained  in  earlier  reports  which  were  delivered  to  the  sponsor  on  three 
occasions:  at  the  Preliminary  Design  Review  (PDR)  on  15  October  1991  (Ref.  1);  at 
an  intermediate  briefing  on  9  January  1992  (Ref.  2),  and  at  the  Critical  Design  Review 
(CDR)  on  14  February  1992  (Ref.  3).  In  the  present  report,  we  discuss  the  methods 
used  in  the  underlying  analysis  and  recapitulate  the  key  results  and  recommendations 
obtained. 

Section  2  discusses  the  goals  and  the  approach  recommended  by  EOS,  Inc.  for 
fielding  data  collection  experiments  in  support  of  laser  Doppler  imaging.  As  detailed  in 
that  Section,  a  two-stage/four-experiment  approach  is  recommended  for  achieving 
these  goals.  Considerations  of  optical  design  and  implementation  are  treated  in 
Section  3.  The  electronic  and  data-col lection  aspects  of  the  experiment  are  treated  in 
Section  4.  Section  5  deals  with  calibration  issues  and  Section  6  discusses  the  Texture 
Imaging  technique.  Section  7  contains  our  summary  and  recommendations. 

1.1  Principal  Accomplishments 

The  initial  period  of  the  program  was  concerned  with  assessing  the  ISTEF  site 
facilities.  These  included  the  seeded,  doubled  Nd:YAG  laser  which  was  physically 
located  at  CREOL,  and  the  Contraves  0.5-m  telescope,  located  at  TTC.  These  two 
major  pieces  of  equipment  are  vital  for  the  experimental  program  leading  to  proof-of- 
coricept  demonstrations  of  autodyne,  Doppler  and  texture  image  tracking  techniques.  ~ 

During  the  course  of  characterizing  this  equipment  and  determining  what  was 
necessary  for  integration  at  ISTEF,  a  sequential  experiment  plan  was  developed.  This 
four  stage  program  provides  an  early  opportunity  to  collect  active  target  returns  to 
verify  photometric  calculations  and  equipment  performance.  The  next  stage  would 
produce  calibrated  measurements  of  target  reflectivity  and  demonstrate  texture 
imaging.  The  third  stage  would  incorporate  interferometric  measurement  of  Doppler 
signatures  to  evaluate  autodyne  potential  and  the  final  stage  would  attempt  closed- 
loop  tracking.  The  plan  allows  for  early  active  data  collection,  with  the  less  complex 
early  stages  designed  for  maximum  likelihood  of  success;  the  plan  should  integrate  - 
well  with  the  development  of  advanced  instrumentation  necessary  for  Doppler  data 
collection,  which  is  necessarily  more  cojnplex. 
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The  work  has  emphasized  the  mechanical,  optical  and  electronic  design  of  the 
apparatus  to  date.  The  Contraves  mount  was  fitted  with  new  coude  mirrors  for  the 
YAG  laser;  the  beam  expander  was  removed  and  a  new  window  installed;  and  a  new 
lens  designed  and  installed  in  the  receiver  to  allow  the  use  of  an  existing  mechanism 
for  focusing.  A  new  back  plane  for  the  Contraves  has  been  designed  to 
accommodate  sensors  for  0.53  pm  and  1.06  pm,  CCD  imaging  and  a  fiber  optic  feed 
to  the  interferometer.  This  back  plane  is  currently  under  construction. 

To  accelerate  the  implementation  of  the  interferometer,  a  fixed,  ruggedized 
version  was  designed  and  ordered.  Field  lenses  for  the  detectors  and  fiber  optics 
were  designed  and  tested  and  the  fibers  themselves  procured.  Isolation  filters  for  the 
CREOL  laser  power  monitors  were  designed,  procured  and  delivered  to  CREOL.  In 
turn,  CREOL  delivered  detectors  for  the  laser  returns  to  EOS.  These  were  tested  and 
were  found  to  require  mechanical  and  electrical  modification  for  compatibility.  Beam 
splitters  and  mounts  to  separate  the  received  signals  by  wavelength  and  function  are 
on  hand. 

Much  of  the  data  recording  capability  resides  at  ISTEF  and  will  be  put  into 
operation  during  the  preliminary  laser  measurement  phase.  EOS  has  developed  the 
software  to  use  a  LeCroy  multichannel  digitizer  to  record  the  detector  outputs. 

Designs  for  synchronization  of  cameras  and  laser  have  been  developed  and  are 
awaiting  implementation. 


2.0  EXPERIMENT  GOALS  AND  APPROACH 


The  objectives  ot  the  experiments  discussed  in  this  report  are  as  follows: 

•  To  obtain  Doppler  data  on  the  plumes  of  in-flight  rockets  by  laser 
interferometric  means-the  first  time  this  will  be  done. 

•  To  obtain  calibrated  data  on  the  plume  reflectance  during  rocket 
flight— also  the  first  time  this  will  have  been  accomplished. 

•  To  demonstrate  the  feasibility  of  the  texture  imaging  concept-another 
first. 

•  To  identify  the  trailing  edge  of  the  rocket  in  an  image  containing 
elements  both  of  the  plume  and  of  the  hard  body.  This  would  be  done 
either  by  means  of  texture-imaging  or  by  combining  the  angle-angle 
image  information  with  Doppler  discrimination  data. 

2.1  Direct-Detection  Doppler  Interferometry  and  Texture  Imaging  Concepts 

The  key  concept  which  these  experiments  seek  to  validate  in  the  field  is  that 
one  can  obtain  supplementary  information  useful  for  tracking  by  discriminating  laser 
returns  scattered  by  the  rocket  plume  from  returns  scattered  by  the  hard  body  through 
the  use  of  suitable  apparatus  added  on  to  a  conventional  telescope.  In  particular,  one 
can  employ  discrimination  based  on  the  Doppler  spectrum.  Alternatively  (or  addition¬ 
ally),  one  may  discriminate  between  the  image  textures  characteristic  of  rockets  and  of 
their  plumes,  as  described  below. 

The  first  method-Doppler  spectra  via  direct  detection-is  illustrated  in  the  upper 
part  of  Figure  1 .  Consider  a  generic  rocket,  traveling  to  the  left,  illuminated  by  a  laser 
spot  whose  diameter  is  comparable  to  that  of  the  rocket.  If  the  spot  reflects  solely 
from  the  rocket  body  (spot  #1  in  the  figure),  a  narrow  Doppler-shifted  return  is  ex¬ 
pected,  as  schematically  illustrated  by  the  inset  figure  at  the  top  left.  If  the  laser  spot 
is  intercepted  by  the  plume  alone  (spot  #3  in  the  figure),  the  Doppler  return  will  be 
centered  at  a  different  frequency,  because  the  mean  of  the  component  of  the  line-of- 
sight  velocity  will  differ  from  that  of  the  rocket,  and  the  spectrum  will  be  spread  out 
because  of  variations  in  the  line-of-sight  velocity  components  of  different  particles  in 
the  illuminated  portion  of  the  plume.  This  characteristic  of  the  plume  Doppler  spec¬ 
trum  is  illustrated  in  the  inset  figure  at  the  top  right.  The  interface  between  the  hard 
body  and  the  plume  will  be  contained  in  the  scatter  from  a  spot  such  as  that  labeled 
#2  in  the  figure.  The  interface  will  not  be  resolvable  from  the  image  alone,  but  can  be 
located  through  the  use  of  the  auxiliary  information  which  is  sought  in  the  experiments 
under  discussion. 

The  conventional  way  to  obtain  the  Doppler  information  from  laser  scatter  return 
signals  would  be  to  use  heterodyne  interferometry,  a  method  best  suited  to  infrared 
wavelengths.  Coherent  detection  methods  are  more  difficult  to  implement  at  visible 
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wavelengths,  whereas  rocket  bodies  do  not  reflect  very  well  in  the  infrared.  The  laser 
Doppler  techniques  of  interest  to  us  employ  direct  detection  instead.  In  the  lower  por¬ 
tion  of  Figure  1 ,  two  such  methods  are  shown.  The  first  is  the  autodyne  method 
shown  on  the  left,  which  yields  relative  Doppler  shifts,  i.e.,  there  is  no  indication  of  the 
sign  of  the  shift.  The  second  technique  is  the  interferometric  method,  which  is  easier 
to  implement  at  the  present  state  of  technology.  In  the  interferometric  approach,  a 
multibeam  Fizeau  interferometer  (MBFI)  is  employed  along  with  anamorphic  optics  to 
generate  the  Doppler  signatures,  and  one  can  distinguish  positive  from  negative 
Doppler  shifts  within  the  free  spectral  range  (FSR)  of  the  instrument.  For  a  typical 
design  considered,  FSR  =  20  to  30  GHz. 

The  other  novel  technique  whose  requirements  we  analyze  is  the  image  texture 
method,  illustrated  schematically  in  Figure  2.  In  the  figure,  we  consider  a  detector  in 
the  form  of  a  linear  imaging  array  of  diffraction-limited  pixels.  Let  the  characteristic 
Doppler  spread  of  the  plume  be  denoted  by  Avp  and  denote  the  Doppler  width  associ¬ 
ated  with  the  hard  body  by  Av  (which  includes  the  finite  coherence  of  the  laser,  the 
effective  resolution  of  the  interferometer,  etc.).  Then,  by  employing  an  integration  time 
T2  such  that  1/Avp  <  T2  <  1/Av  (e.g.,  T2  on  the  order  of  100  ns),  one  can  distinguish 
between  hard  body  scattering  and  plume  scattering  by  the  differently  fluctuating 
natures  of  the  two  types  of  returns:  the  plume  image  will  tend  to  have  a  smoother 
distribution  of  intensity  vs  pixel  position  than  will  the  highly  speckled  image  expected 
from  the  hard  body. 

2.2  Target  Kinematics  at  ISTEF 

The  ISTEF  site  in  Florida  is  uniquely  well  suited  for  tracking  rocket  bodies  and 
their  plumes  during  the  period  immediately  following  ignition  and  liftoff.  For  purposes 
of  designing  an  experiment  to  test  laser  radar  tracking  concepts  at  ISTEF,  we  consider 
the  Delta-2  rocket,  whose  essential  mechanical  features  are  indicated  in  Figure  3. 

The  essentially  cylindrical  rocket  body  is  30  to  40  meters  long;  it  is  powered  by  rocket 
motors  whose  outer  periphery  is  contained  in  a  circle  4.5  m  in  diameter.  The  plume 
properties  during  the  lift-off  episode  visible  from  ISTEF  are  governed  by  the  propellant 
characteristics  listed  in  the  figure;  scattering  contributions  are  to  be  expected  from  the 
exhaust  plumes  of  both  types  of  rocket,  i.e.,  liquid-  (LOX/RP-1)  and  solid-fueled  (14%' 
Al). 


Based  on  consideration  of  the  rocket  kinematics,  one  can  estimate  the  quantity 
of  data  to  be  expected  during  the  course  of  a  lift-off  experiment.  For  the  Delta  rocket, 
with  a  thrust  T  =  860  klbs  and  weight  W  =  484  klbs,  the  thrust  length  will  be  approxi¬ 
mately  40  m  and  the  target  acceleration  will  be  a  =  (T/W  -1)g  =  7.78  ms’2,  so  that  the 
target  velocity  vT  (m/s)  will  increase  with  altitude  h  (m)  as  vT  =  (2ha),/z .  We  consider 
an  experiment  in  which  the  laser  illumination  forms  a  spot  aimed  at  a  fixed  altitude 
and  the  target  flies  through  the  beam.  The  projected  laser  spot  will  have  a  beam 
diameter  DB  v*n);  in  order  to  "paint"  the  entire  rocket  body  with  such  spots  at  a  pulse 
repetition  frequency  PRF  (Hz)  with  an  overlap  factor  OL  between  adjacent  spots,  DB 
must  satisfy  the  following  relationship:  .. 
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Db  =  (2ha)1/2  (1+OL)/PRF. 


Taking  PRF  =  10  Hz,  the  spot  diameter  projected  by  the  fixed  laser  beam  should  lie 
between  the  two  limiting  curves  indicated  on  Figure  4  in  order  to  cover  the  target 
more-or-less  uniformly  during  the  data  collection  episode.  The  two  limiting  curves 
assume  either  (1)  no  overlap  between  successive  beam  positions,  or  (2)  50%  overlap 
between  successive  spots.  As  noted  on  the  figure,  the  spot  size  required  (from  kine¬ 
matic  considerations  alone)  becomes  larger  than  the  width  of  the  booster  for  altitudes 
above  60  m  or  so. 

The  length  of  time  that  the  target  remains  in  the  field  of  view  in  such  an  experi¬ 
ment  can  now  be  calculated.  Hence,  as  Figure  5  indicates,  the  number  of  data 
frames  associated  with  the  rocket  body  as  it  passes  through  the  fixed  laser  beam  will 
be  on  the  order  of  10  to  30  (depending  on  beam  altitude),  if  one  assumes  that  the 
data  rate  is  governed  by  the  10-Hz  PRF.  Thus,  an  initial  experiment  predicated  upon 
having  the  rocket  pass  through  a-fixed  laser  beam  should  be  capable  of  yielding 
statistically  significant  data,  provided  that  the  rocket  does  pass  through  the  beam. 

2.3  Step-wise  Approach  to  Data  Collection 

From  the  results  shown  in  Figures  4  and  5,  the  constraints  imposed  by 
requiring  the  laser  to  have  a  fixed  pointing  direction  become  clear.  To  obtain  data  on 
rocket  bodies  and  plumes  at  altitudes  above  60  m  or  so,  it  will  be  necessary  to  go 
beyond  such  a  first-phase  experiment.  Of  course,  once  tracking  of  the  target  during 
the  experiment  is  required,  more  complex  hardware  will  be  needed. 

EOS,  Inc.  has  given  serious  consideration  to  the  nature  of  the  information  which 
it  is  reasonable  to  expect  can  be  obtained  with  and  without  tracking.  Figure  6-- 
originally  presented  at  the  PDR  (Ref.  1 )  and  repeated  at  the  CDR-summarizes  the 
two-stage,  four-experiment  step-wise  approach  for  data  collection  which  we  have 
developed.  The  two  experiments  during  Stage  1  are  limited  to  the  lift-off  episode  and 
do  not  obtain  Doppler  signature  data  and  do  not  require  tracking  the  target.  Experi¬ 
ments  3  and  4,  both  aimed  at  obtaining  Doppler  rignature  data,  would  have  to  take 
place  during  Stage  2,  since  the  target  would  need  to  be  tracked.  In  Experiment  3, 
manual  tracking  would  suffice,  but  automatic  means  of  tracking  would  be  highly 
desirable  for  Experiment  4. 

2.3.1  Experiment  1:  No  Tracking,  2-m  Spot 

The  data  collection  plan  for  Experiment  1  calls  for  obtaining  calibrated  returns 
from  the  laser-illuminated  hard  body  and  target  at  two  wavelengths,  one  in  the  visible 
(doubled  Nd:YAG  at  0.532  pm)  and  one  in  the  near-infrared  (Nd:YAG  at  1.063  pm). 

The  experiment  plan  also  calls  for  demonstrating  the  ability  to  perform  spatially 
resolved  texture  imaging  at  the  shorter  wavelength  (in  the  laboratory,  in  preparation  for- 
Experiment  2). 
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To  implement  Experiment  1  (Figure  7),  one  would  employ  a  laser  illuminator 
pointing  at  a  fixed  location,  and  no  tracking  would  be  required  (although  manual 
tracking  would  be  valuable,  of  course,  as  would  real-time  range  data).  The  optical 
design  would  yield  a  spot  diameter  of  approximately  2  m  at  the  target,  without  the 
need  for  beam-forming  optics.  The  target  vehicle  would  move  through  the  fixed  beam 
while  the  receiver  gathers  the  scattered  return.  The  receiver  design  calls  for  a  50-cm 
collecting  aperture,  with  the  received  signals  split  into  three  paths  (Figure  8  (a)).  The 
first  path  leads  to  an  infrared  "photon  bucket"  detector  such  as  an  avalanche 
photodiode  (APD).  The  second  path  leads  to  an  intensified  charge-coupled  detector 
(ICCD)  camera,  gated  synchronously  to  the  laser  pulse  to  achieve  background 
rejection,  especially  self-luminous  and  ambient-scattered  radiation  from  the  plume. 

The  third  path  leads  to  a  "photon  bucket"  (e.g.,  an  APD)  for  the  visible  channel  at 
0.532  pm. 

2.3.2  Experiment  2:  No  Tracking,  5-m  Spot 

The  data  collection  plan  for  Experiment  2  is  designed  to  obtain  laser  returns 
from  complex  targets,  including  glints,  etc.  Emphasis  will  be  placed  on  expanding  the 
calibrated  target  cross-section  data  base  at  the  same  wavelengths  as  in  Experiment  1 . 
Another  major  goal  is  to  obtain  spatially  resolved  texture  images  of  the  target  at  0.532 
pm.  (As  an  alternative,  should  it  be  judged  that  sufficient  target  cross-section  data  is 
available  following  Experiment  1 ,  Doppler  interferometric  measurements  could  be 
commenced  during  Experiment  2.) 

To  implement  Experiment  2  (Figure  9),  modifications  to  the  optics  would  be 
desirable,  so  as  to  increase  the  beam  size  while  matching  the  detector  field  of  view  to 
the  beam  (Figure  8(c)).  Beam-forming  optics  would  allow  the  beam  size  to  increase  to 
5  m  at  in  the  infrared.  (Should  beam  forming  optics  not  become  available  in  a  timely 
manner,  beam  expanders  would  be  omitted,  and  the  MBFI  would  be  used  with  the 
same  green  laser  line  as  in  Experiment  1 .)  As  was  the  case  for  Experiment  1 , 
Experiment  2  could  be  performed  with  the  laser  beam  fixed  in  position,  but  it  would  be 
highly  desirable  to  be  able  to  track  the  target  manually.  For  this  purpose,  a  pre¬ 
programmed  trajectory  would  be  required,  along  with  real-time  range  updates. 

Finally,  as  indicated  in  Figure  8(b),  should  the  signal  levels  expected  at  the 
photodetectors  be  sufficiently  high,  one  would  replace  the  ADP  at  0.532  pm  (in  the 
third  path  in  Figure  8  (a))  by  an  optical  fiber  connecting  to  the  MBFI  followed  by  a 
gated  camera. 

2.3.3  Experiment  3 

Experiment  3  calls  for  obtaining  Doppler  data,  from  both  the  hard  body  and  the 
plume,  to  evaluate  the  potential  of  autodyne  tracking  techniques.  To  implement 
Experiment  3,  the  green  APD  would  be  replaced  by  an  optical  fiber  receiver  with  a 
field  lens  whose  numerical  aperture  is  matched  to  th.  of  the  fiber.  The  other  end  of 
the  fiber  will  be  located  in  the  laboratory  and  its  output  will  feed  the  MBFI. 
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he  results  of  Experiments  1  and  2  will  be  used  to  determine  the  signal  levels 
that  ca,  be  expected,  so  that  the  optimum  split  between  the  CCD  camera  and  the 
MBFI  can  be  selected  and  the  appropriate  beam  splitter  installed  in  the  receiver. 
Beam  splitters  with  a  range  of  splitting  ratios  from  10  percent  to  50  percent  are  on 
hand. 


As  described  in  Section  3.1,  there  is  a  set  of  optimum  parameters  for  the 
interferometer.  Inasmuch  as  the  interferometer  is  located  near  the  exit  pupil  of  the 
receiver  telescope,  the  telescope  field-of-view  will  be  restricted.  For  full  utilization  of 
all  green  photons  in  the  interferometer,  the  optimum  field-of-view  would  be  0.2  milli- 
radian,  which  corresponds  to  a  2-meter  spot  at  a  target  10-km  distant.  This  would 
necessitate  adjustment  of  the  beamforming  optics  from  the  5-meter  spot  in  Experiment 
2.  Depending  on  the  results  of  the  laser  characterization  tests  and  the  consequent 
beam  expander  design,  it  is  planned  to  maintain  the  IR  beam  as  in  Experiment  2. 

2.3.4  Experiment  4 

At  present,  ISTEF  does  not  employ  closed-loop  target  tracking.  For  this 
reason,  during  Experiments  1, 2  and  3,  the  illumination  beam  is  initially  fixed  and 
directed  above  the  target,  pre-launch,  and  the  target  is  allowed  to  move  through  the 
beam.  Tracking  will  be  done  manually  after  the  missile  and  plume  have  moved 
through  the  beam. 

It  is  expected  that,  by  the  time  Experiment  4  is  implemented,  some  form  of 
closed-loop  tracking  will  be  exercised.  It  is  intended  to  build  into  the  tracking  routine  a 
programmed  variable  bias  so  that  the  illuminating  beam  will  be  "painted”  over  the 
target.  The  specific  details  depend  on  what  track  algorithm  is  employed.  This  will  be 
developed  during  the  course  of  closed-loop  track  implementation.  If  the  closed-loop 
track  can  follow  the  target  with  a  substantial  increase  in  the  target  range,  it  may  be 
necessary  to  exercise  other  variables.  These  include:  a  sliding  range  gate; 
programmed  beamforming  optics  to  maintain  beam  size;  and  programmed  focus 
control.  Inputs  for  controlling  these  variables  can  be  derived  from  predictive  range 
estimation  and  updated  by  real-time  range  measurement  using  the  laser  pulse  and  a 
timer.  It  is  intended  to  investigate  this  capability  during  the  manual  track  phase  of 
Experiments  1  and  2. 
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3.0  OPTICAL  DESIGN  AND  IMPLEMENTATION 


The  optical  design  for  Experiments  1  and  2  was  performed  in  sufficient  detail  to 
draw  up  a  hardware  list  for  presentation  at  the  CDR  (Figure  10).  Several  key 
components  were  ordered  and  received,  as  listed  in  Figure  11.  As  received,  the 
avalanche  photodiodes  (APD’s)  were  not  suitable  and  they  have  had  to  be  modified. 
The  optical  parameters  of  the  field  lenses  for  the  various  detectors  are  listed  in  Figure 
12.  The  beam  splitters,  which  were  designed  for  a  22.5-deg  incidence  angle,  have 
been  tested.  The  measured  values  for  the  splitting  ratios  are  compared  with  the 
design  values  in  Figure  13. 

3.1  Design  of  the  Interferometer  and  Transfer  Optics 

The  MBFI  design  calls  for  an  optically  contacted  air-spaced  wedge  etalon  to 
ensure  reliability  (Figure  14).  RFQ’s  were  issued  to  four  potential  vendors,  eliciting 
two  responses  (from  Continental  and  from  TecOptics).  Continental’s  bid  was 
accepted.  Their  12-week  delivery  time  (measured  from  18  February)  should  allow 
testing,  assembly  and  shipment  in  time  for  availability  at  ISTEF  in  July. 

A  key  parameter  in  designing  the  transfer  optics  linking  the  MBFI  to  the  receiver 
aperture  is  the  field  of  view  required  for  experiments  at  ISTEF.  As  indicated  in  Figure 
1 5,  the  maximum  angular  deviation  from  the  nominal  line  of  sight  is  magnified  by  the 
aperture  ratio  in  a  matched-f/#  design.  The  reason  for  concern  is  indicated  by 
consideration  of  the  data  shown  in  Figures  16  and  17.  The  calculations  summarized 
in  these  figures  assumed  the  nominal  parameters  for  the  MBFI  shown  on  the  upper 
left:  Free  spectral  range  FSR  =  15  GHz  at  wavelength  X=  0.532  pm;  equal  93% 
reflectances  (essentially  lossless)  at  the  two  surfaces  which  are  inclined  at  a  wedge 
angle  a  =  1 6  prad. 

The  abscissa  in  Figure  16  gives  the  nominal  incidence  angle  0O  with  respect  to 
the  normal  to  the  first  reflecting  surface;1  the  calculations  assumed  a  field-of-view 
subtending  a  cone  of  half-angle  0S=  1  mrad  (i.e.,  FOV  =  2  mrad).  The  triangle-bearing 
dashed  curve  in  Figure  1 6  plots  the  variation  in  peak  transmission  with  incidence 
angle  (right-hand  scale),  assuming  that  fringe  intensity  contributions  associated  with 
collimated  light  originating  from  different  directions  add  incoherently;  a  9-point 
quadrature  rule  was  used.  When  the  line  of  sight  deviates  by  more  than  =1  mrad  from 
the  optimal  value  (0o=-O.5  mrad),  the  peak  intensity  drops  two-fold. 

The  solid  curve  in  Figure  16  indicates  the  location  of  the  average-intensity  peak, 
while  the  short-dashed  curves  indicate  the  locations  of  the  half-intensity  points  (relative 


1  A  negative  sign  indicates  incidence  from  the  side  of  the  normal  opposite  the  apex  of  the 
wedge.  The  calculation  was  done  for  fringes  observed  just  past  the  second  reflecting 
surface.  For  fringe  observation  in  other  planes,  the  optimal  incidence  angle  changes,  but 
the  curves  in  Figure  16  retain  their  shape. 
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to  the  peak  intensity  for  a  given  0O).  Fringe  location  is  expressed  in  milli-orders  (m- 
order)  on  the  left-hand  ordinate  scale,  i.e.,  in  units  of  FSR/1000;2  since  FSR  =  15 
GHz,  a  shift  of  10  m-order  corresponds  to  a  frequency  shift  of  150  MHz.  The  curves 
in  Figure  16  indicate  that  a  field  lens  designed  for  a  1-mrad  field  of  view  should  keep 
the  FWHM  below  *300  MHz,  and  that  a  well-aligned  MBFI  can  locate  peaks  to  an 
accuracy  <100  MHz. 

The  choice  of  mirror  reflectances  for  the  interferometer  plates  and  the  tolerance 
to  be  placed  on  the  field  of  view  have  been  determined  by  means  of  the  calculations 
such  as  the  one  summarized  in  Figure  17.  For  this  calculation,  the  reflectances  of  the 
plates  are  94%  and  the  spacing  between  the  plates  has  been  increased  so  that  FSR 
is  now  12  GHz  (as  indicated  in  the  figure  heading).  The  solid  curve  in  Figure  17 
indicates  the  dependence  of  the  FWHM  (full  width  between  half-transmission  points 
near  the  peak  intensity-left-hand  ordinate  scale)  upon  the  field  of  view,  i.e.,  as  a 
function  of  the  half-angle  of  the  cone  defining  the  field  of  view  (FOV  =  20s).  The 
dashed  curve  plots  the  field-of-view  dependence  of  the  maximum  transmission  of  the 
interferometer,  averaged  over  the  field  of  view  (to  be  read  from  the  right-hand  scale). 
With  this  design,  a  maximum  (averaged)  transmission  of  80%  or  better  requires  a 
smaller  field  of  view  (20s  =  1.5  mrad  instead  of  2  mrad).  The  associated  FWHM  would 
then  be  greater  (=500  MHz)  despite  the  decreased  free  spectral  range. 

Figure  18  shows  the  top-level  considerations  in  designing  transfer  optics 
employing  a  fiber  optic  light  guide.  The  focal  length  is  related  to  the  fiber  diameter  by 
the  considerations  shown  in  Figure  19.  In  our  design,  we  first  characterized  available 
fiber  optics  by  their  numerical  apertures  at  the  input  and  output  ends.  Then  we 
selected  a  suitable  combination  of  fiber  diameter  and  numerical  aperture,  selected  a 
field  lens  to  go  with  it.  The  next  step  was  to  design  the  collimator. 

3.2  Beam  Expander  Design 

As  indicated  above  (Sec.  2.3.2),  it  is  desirable  for  a  properly  designed  beam 
expander  to  be  available  in  time  for  Experiment  2.  This  experiment  is  designed  to 
acquire  signature  data  at  two  wavelengths,  one  in  the  infrared-with  a  5-m  beam-and 
the  other  in  the  visible-with  a  smaller  beam  diameter  needed  for  the  Doppler 
measurements.  A  beam  expander  is  desirable  to  maintain  control  of  the  optical  paths 
independently  at  the  two  wavelengths.  The  optical  design  must  ensure  proper  optical 
correction  at  each  wavelength,  while  maintaining  collinearity  of  the  two  beams  (Figure 
20).  The  philosophy  underlying  the  design  of  the  beam  expander  is  summarized  in 
Figure  21 . 


2  The  "zero"  reference  is  the  expected  position  of  the  peak  at  normal  incidence;  it  is  the 
peak  position  only  in  the  (Fabry-Perot)  limit  that  the  wedge  angle  tends  to  zero. 
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3.3  Hardware  Summary 


The  following  salient  points  were  made  at  the  Critical  Design  Review,  14 
February  1992: 

•  The  design  and  testing  of  hardware  for  Experiment  1  is  on  schedule. 

•  The  optical  design  for  the  autodyne  and  Doppler  interferometric 
measurements  achieves  an  instrument  performance  which  is  within  a 
factor  of  2  of  the  theoretical  optimum  based  on  first  principles. 

•  Resolution  of  remaining  design  uncertainties  are  not  expected  to 
cause  delay  in  data  collection. 

•  Definitive  measurements  of  the  laser  energy  distribution  are  needed  in 
order  to  ensure  photometric  integrity  of  the  Doppler  data. 

•  Authorization  from  ONR  to  proceed  with  procurement  was  sought  and 
obtained  at  the  CDR. 
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4.0  ELECTRONIC  DAT A-COLLECTION  INTERFACE 


Figure  22  constitutes  a  block  diagram  of  the  control  and  signal  circuitry 
associated  with  the  telescope  and  the  apparatus  to  be  mounted  on  it  for  Experiment  1 . 

Since  several  organizations  shall  be  participating  in  the  experiment,  it  is 
essential  that  the  responsibilities  of  each  be  spelled  out.  The  break-out  of 
responsibilities  proposed  at  the  CDR  is  as  follows: 

(1)  EOS,  Inc.  would  be  responsible  for  the  following: 

•  Data  acquisition  and  control  computer  and  its  associated  hardware 

•  CAMAC 

•  Avalanche  photodiodes  and  associated  DC  power  supplies  and 
amplifiers. 

•  Signal  compatibility-from  the  laser  trigger  to  the  CAMAC,  from  the 
time  gate  generator  to  the  Xybion  camera,  and  from  the  laser  power 
supply  to  the  CAMAC. 

•  Subject  to  Air  Force  approval,  EOS  would  employ  the  Cohu  camera 
(owned  by  the  USAF)  for  monitoring  the  beam  quality. 

(2)  NRC/ISTEF  would  be  responsible  for  the  following: 

•  Provide  the  Xybion  and  tracking  cameras.  Using  synchronization 
signals  derived  from  the  Xybion  camera,  gen  lock  the  tracking  and  beam 
monitoring  camera. 

•  Synchronize  the  laser  triggering  to  the  Xybion  camera: 

•  Obtain  a  frame  counter  (a  30  Hz  signal  indicating  the 
start  of  a  frame)  from  the  Xybion  camera  or  any  other 
device  that  uses  the  camera’s  sync  signals  as  input. 

•  Divide  the  frame  counter  rate  by  3  to  obtain  the  desired 
laser  pulse  repetition  frequency  of  1 0  Hz. 

•  Send  the  10-Hz  signal  to  the  laser  to  initiate  charging  of 
the  laser  lamps. 

•  Sense  the  laser’s  end-of-charge  signal. 

•  Send  the  "fire"  command. 

The  required  pulse  durations  and  amplitudes  are  listed  in  the  laser 
manual  or  can  be  obtained  directly  from  the  manufacturer  (Continuum). 

•  Make  available  the  Tektronix  Model  602  digitizing  signal  analyzer. 

•  Record  all  video  signals  with  IRIG  B 
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•  Supply  the  IRIG  B  signal  to  EOS 

•  Supply  cabling  between  remote  locations  including: 

•  2  RG-58  cables  with  BNC’s  at  least  200  ft  long 

•  Cables  from  the  sensors  to  the  recording  locations. 

•  Sync  cable  to  the  laser. 

•  Provide  means  (via  relay  or  otherwise)  of  starting  the  laser 
immediately  prior  to  rocket  ignition. 

•  Supply  video  monitors  where  needed. 

•  Supply  intercoms  and/or  other  communication  equipment  as  needed. 

•  Supply  100  VAC  power  for  EOS  equipment;  two  15-amp  outlets 
should  suffice. 

(3)  CREOL  would  be  responsible  for  the  following  hardware: 

•  Nd:YAG  laser  (as  specified  elsewhere). 

•  Laser  power  monitors,  incorporating  filters  supplied  by  EOS. 
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5.0  CALIBRATION  ISSUES 

Issues  relating  to  radiometric  calibration  for  field  experiments  at  ISTEF  were 
treated  in  Refs.  1  and  2  and  are  summarized  here.  Briefly,  calibration  is  essential  (1 ) 
to  ensure  the  colinearity  of  the  optical  systems  for  the  visible  (0.53  pm)  and  infrared 
(1 .06  pm)  beams,  as  well  as  (2)  to  check  for  significant  departures  of  the  illumination 
patterns  from  their  nominal  far-field  shapes. 

The  method  suggested  for  ensuring  colinearity  is  to  employ  ancillary  sensors  in 
both  the  visible  and  infrared  regions  which  observe  a  target  calibration  board 
containing  reference  beacons.  Each  camera  would  be  aligned  with  the  center  beacon; 
spot  locations  would  be  measured  with  respect  to  a  reticle.  Additional  beacons  could 
be  used  to  estimate  r0  so  as  to  quantify  optically  significant  atmospheric  turbulence 
conditions.  If  done  carefully,  the  calibration  procedure  also  may  be  able  to  estimate 
the  atmospheric  transmission  along  the  path.  (See  Figure  23,  from  Ref.  1 ). 

The  fields  of  view  of  the  various  optical  devices  considered  for  use  in 
Experiment  1  are  summarized  in  Figure  24  (from  Ref.  2);  the  fields  of  view  of  the 
Contraves  and  the  Xybion  camera  are  shown  as  centered  on  a  Delta-2  rocket  body, 
viewed  from  a  range  of  10  km.  Full  field  for  the  Xybion  corresponds  to  16  mm  along 
the  diagonal;  the  corresponding  diagonal  for  the  ITT  camera  is  18  mm.  The  dashed- 
line  rectangles  indicate  the  maximum  fields  (along  the  vertical  and  horizontal  direc¬ 
tions)  when  imaging  the  target  with  one  of  three  lenses  on  the  ITT  camera:  10"  f/12 
(4.7  mr  x  3.5  mr),  12"  f/12  (3.9  mr  x  3.0  mr),  and  12"  f/15  (3.1  mr  x  2.4  mr). 

Figure  25  (from  Ref.  2)  describes  one  method  of  system  calibration  based  on 
the  imaging  of  a  reference  beacon.  Conceptually,  one  would  locate  the  center  of  the 
field  of  view  relative  to  the  beacon  image,  as  indicated  in  the  figure  for  the  example  of 
the  ITT  camera  employing  the  12"  f/15  lens. 

The  target  reflectance  in  the  green  region  should  be  on  the  order  of  0.3/sr, 
measured  without  filters,  at  night.  As  indicated  in  Figure  26  (from  Ref.  1),  one  could 
place  a  painted  target  board  on  the  rocket  gantry  and  calibrate  immediately  prior  to  lift¬ 
off.  The  placement  of  a  target  board  on  the  rocket  gantry  may  not  be  practical.  An 
alternative  site,  a  lighthouse  at  nearly  the  same  range,  with  a  similar  azimuth  angle, 
has  been  identified.  This  site  might  be  suitable  for  calibration  once  the  reflectance 
characteristics  have  been  measured. 
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6.0  TEXTURE  IMAGING  PERFORMANCE  ANALYSIS 


An  introduction  to  the  image  texture  method  was  presented  in  Section  2  (see 
Figure  2).  The  purpose  of  a  field  experiment  for  texture  imaging  is  to  demonstrate  the 
applicability  of  the  technique  for  obtaining  bounds  on  the  plume-induced  Doppler  shifts 
and  for  locating  the  trailing  edge  of  the  rocket.  EOS,  Inc.  has  analyzed  the  feasibility 
of  performing  a  demonstration  experiment  on  texture  imaging  at  ISTEF.  The  analysis 
consists  of  the  following:  (a)  signal-to-noise  (radiometric)  calculations,  (b)  estimation  of 
the  effects  of  clear  air  turbulence  on  the  image  texture  (leading  to  an  expression  for 
speckle  contrast),  and  (c)  performance  predictions  for  the  conditions  likely  to  apply  at 
ISTEF. 

6.1  Radiometric  Analysis 

In  the  measurement  of  target  reflectance,  a  key  quantity  is  the  signal-to-noise 
ratio  of  the  measurement.  We  may  model  the  total  number  ntolal  of  photocounts 
received  in  terms  of  an  overall  system  constant  Csysl ,  along  with  the  beam-target 
overlap  function  \\f  (which  includes  the  effects  of  jitter),  the  two-way  atmospheric 
attenuation  xa ,  the  effective  target  reflectance  p  and  an  additive  term  nb 
(representing  background  and  electronic  noise)  as  follows: 

^total  =  ^syst  V  P  + 

Reflectance  estimates  are  based  on  mean  values,  i.e., 

P  =  <  ntotaI  >  /  Csyst  <  \j/  >  <  xa  >. 

However,  since  photocounts  fluctuate,  we  must  also  consider  the  variance  of  the  total 
number  of  photocounts: 


var  (ntoCal)  _  yar(q r)  +  var(xa)  +  i 

l<ntotal>-<nb>l2  <Ta>2  SNRbkgd 


where 


SNRbkffd  =  [<ncoCal>  -  <nb>]2/var(nb ) 


(2) 


Table  I  indicates  the  target  and  system  parameters  employed  for  radiometric 
calculations  by  EOS,  as  well  as  the  parameters  employed  in  independent  calculations 
by  P.  Gatt  of  CREOL  and  by  F.  Tart  of  W.J.  Shafer  &  Associates.  Similarly,  Table  II 
lists  the  assumptions  made  concerning  the  plume  background  radiometry  by  EOS  and 
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by  the  same  two  other  parties.  Justification  for  the  assumptions  made  by  EOS  has 
been  documented  in  Ref.  4.  (See  also  Refs.  5  and  .)  An  Appendix  to  this  report 
contains  the  equations  used  for  the  radiometric  calculations.  Here,  it  suffices  to 
present  the  results  of  these  calculations. 

Figure  27  shows  the  total  number  of  signal  photons  collected  as  a  function  of 
range  (upper  two  curves)  which  are  associated  with  the  plume  and  with  the  hard  body, 
respectively.  (As  noted  in  Table  I,  we  assume  that  the  plume  reflectance  is  twice  that 
of  the  hard  body.)  For  these  calculations,  the  beam  diameter,  which  is  2  m  at  a  range 
of  10  km  (as  in  Experiment  1),  scales  linearly  with  range.  When  the  total  number  of 
photons  collected  is  divided  by  the  number  of  resolution  cells  in  the  image  of  the 
plume,  the  lower  (solid)  curve  is  obtained,  representing  the  mean  number  of 
photocounts  per  resel  in  the  plume,  for  slant  ranges  appropriate  to  launches  likely  to 
be  visible  from  ISTEF  in  Experiment  1;  one  expects  to  have  on  the  order  of  100  to 
1000  signal  photocounts  per  resel  at  all  realistic  ranges. 

When  the  calculation  is  repeated  for  the  5-m  beam  diameter  appropriate  to 
Experiment  2,  Figure  28  is  obtained.  (The  line  types  have  the  same  significance  as  in 
Figure  27.)  At  ranges  less  than  17  km,  one  still  expects  in  excess  of  100  signal 
photocounts  per  resel  on  the  plume. 

Turning  now  to  the  photocount  level  expected  from  the  plume  background,  we 
obtain  the  results  shown  in  Figure  29  for  two  different  assumptions.  In  case  1 ,  the 
receiver  footprint  at  10  km  is  3  m,  as  would  be  appropriate  to  Experiment  1.  Case  2, 
which  assumes  a  6-m  footprint,  corresponds  to  Experiment  2.  Since  the  plume  fills 
the  field  of  view,  and  since  atmospheric  turbulence  is  expected  to  limit  the  achievable 
resolution  for  this  viewing  geometry,  both  cases  produce  the  same  number  of 
background  photocounts  per  turbulence-limited  resolution  cell.  For  this  first-cut 
calculation,  a  constant  value  of  5  cm  is  assumed  for  the  coherence  diameter  r0 , 

independent  of  range,  and  the  area  of  the  associated  resolution  cell  is  *  (—) 2 . 

T  r° 

(Further  discussion  of  atmospheric  turbulence  is  deferred  until  Sec.  6.2.) 

Next,  we  combine  the  results  shown  in  the  preceding  two  figures  to  obtain  the 
"naive"  lidar  signal-to-noise  ratio  for  observing  the  plume  signal  against  the  plume 
background  noise-SNRy^  of  Eq.  (2)  (i.e.,  ignoring  all  other  sources  of  noise  in  Eq. 
(1))--as  shown  in  Figure  30.  Even  for  case  2  (i.e.,  Experiment  2),  the  background- 
noise  limit  to  the  per-resel  signal  to  noise  ratio  exceeds  10:1  at  all  ranges,  and 
exceeds  100:1  at  ranges  below  18  km. 

6.2  Speckle  Contrast  Analysis 

We  showed  many  years  ago  (Refs,  7-8)  that  one  may  employ  a  linear-systems 
approach  to  calculation  of  the  speckle  contrast  of  laser  scatter  from  a  rough  target  in  *" 
the  presence  of  atmospheric  turbulence.  As  indicated  in  Figure  31 ,  we  consider  a 
transmitter  aperture  at  plane  X,  which  Illuminates  a  target  at  plane  X„ .  The  scattered 
radiation  is  observed  by  the  receiver  at  plane  X’,  (co-located  with  the  transmitter) 
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which  images  the  target  onto  the  image  plane  X2 .  At  the  image  plane,  the  complex 
amplitude  associated  with  the  k'h  scatterer,  denoted  by  Uk  (Xj, )  is  calculable  as  the 
triple  convolution  indicated  on  the  figure.  The  notation  is  such  that  B(x0-x1 ) 
represents  the  "raw"  beam  amplitude  at  the  target  (in  the  absence  of  turbulence),  GA 
(•)  represents  a  Green’s  function  associated  with  the  turbulence,  h(»)  denotes  the 
transfer  function  of  the  receiving  telescope,  and  the  target  reflectance  function 
possesses  a  time-dependent  argument  which  associates  the  complex  amplitude  with 
the  klh  scatterer.  The  intensity  measured  by  the  detector  is  calculated  from  the 
squared  modulus  of  Uk ,  summed  over  all  scatterers  in  the  region  of  overlap  between 
the  beam  footprint  and  the  receiver  field  of  view.  We  make  a  delta-correlation 
(localization)  approximation  such  as  is  conventional  for  analysis  of  speckles  in 
scattering  from  rough  surfaces.  The  final  result  of  a  triple-product  covariance  analysis 
is  that  the  speckle  contrast  will  not  be  much  perturbed  by  the  presence  of  turbulence, 
provided  that  the  coherence  diameter  r0  is  smaller  than  the  (turbulence-limited) 
resolution  cell. 

We  follow  the  example  of  David  Fried  (Ref.  9)  in  representing  the  influence  of 
clear-air  turbulence  on  the  optical  system  by  the  coherence  diameter,  ra ,  which  is 
defined  on  the  basis  of  the  wave  structure  function  Dw(r),  where  r  denotes  trans¬ 
verse  separation  across  the  wave-front;  this  function  represents  the  mean-square  devi¬ 
ation  of  the  (complex)  wavefront  from  the  value  it  would  have  in  the  absence  of  turbu¬ 
lence.  In  the  Kolmogorov  inertial  subrange  of  turbulence  (well-developed  turbulence), 
the  wave  structure  function  scales  with  the  coherence  diameter  according  to  the  5/3- 
law: 


DJr)  =  6.44  (i/f#)“ 

where  ra  is  a  function  of  the  wavelength  and  of  the  altitude-  and  range-dependent 
micrometeorology.  For  homogeneous  turbulence,  the  range-  and  wavelength-depen¬ 
dences  can  be  factored  out;  for  more  complicated  situations,  numerical  integration  is 
required.  For  purposes  of  estimating  rol  we  started  with  the  turbulence  profile  sug¬ 
gested  by  Greenwood  (Ref.  10)  (which  was  based  on  data  at  the  ARPA  Maui  Optical 
Station  (AMOS),  Ref.  1 1 )  and  adapted  it  to  the  near-sea-level  site  geometry  at  ISTEF- 
For  a  slant  range  R  =  10  km  between  points  located  at  heights  h  =  5  to  10  m  above 
MSL,  the  model  predicts  the  following  range  of  values  for  ra: 

r0  =  5  to  16  mm  at  0.532  pm 
=  1 2  to  36  mm  at  1 .064  pm 

Such  values  are  consistent  with  data  obtained  near  ground  level  in  carefully  controlled 
field  measurements  conducted  by  RADC  over  grassy,  more-or-less  level  terrain  near 
Verona,  NY  (Ref.  12).  Altitude  dependence  is  assumed  to  follow  a  -1.3  power  law, 
i.e.,  turbulence  strength  scales  according  to  (h+h0)'13,  with  h0=  0.01  m.3 


3  To  be  precise,  we  took  Cn2  =[.0043  exp(-h/4000.)+  22.(h+.01)'13]x10'14  [m'2*]. 
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For  present  purposes,  we  estimate  the  turbulence-limited  effective  resolution  of 
the  receiver  telescope,  Atoa,  from  the  approximate  relation: 


A*.  *  Mr0 


where  X  is  the  wavelength  and  R  is  the  slant  range  (mks  units). 

Figure  32  displays  the  results  of  applying  the  turbulence  model  to  computing 
the  coherence  diameter  and  the  turbulence-limited  receiver  resolution  (both  expressed 
in  mm)  at  X=0.532  pm,  over  horizontal  ranges  <10  km,  between  points  at  altitude  5, 

10  and  30  m  above  sea  level.  It  is  evident  from  the  plots  that  the  coherence  diameter 
is  likely  to  be  significantly  smaller  than  A^  at  all  realistic  ranges  available  at  ISTEF. 

It  is  important  to  take  the  turbulence  levels  properly  into  account;  if  one  were  to 
employ  uncharacteristically  small  values  for  the  Kolmogorov  structure  parameter  near 
ground  level,  one  would  obtain  results  such  as  illustrated  in  Figure  33,  where  the 
same  power  scaling  with  (h+.01[mj)'13  is  assumed,  but  C„2  at  h=10  m  has  the  smaller 
value  indicated  on  the  figure. 

At  X  =  1 .064  pm,  the  expected  coherence  diameter  and  resolution  cell  sizes 
are  as  shown  in  Figure  34--for  same  the  realistic  turbulence  model  as  in  Figure  32. 
Although  at  this  longer  wavelength  the  coherence  diameter  is  greater  and  the  turbu¬ 
lence-limited  resolution  cell  is  smaller,  it  is  still  expected  that  r0  <  Ares  at  ranges 
exceeding  =5  km. 

A  related  issue  in  an  image  texture  experiment  is  the  amount  of  contrast  modu¬ 
lation  that  one  expects  on  account  of  the  beam  steering  induced  by  atmospheric  tur¬ 
bulence.  An  approximate  analytical  model  for  the  contrast  modulation  induced  by 
such  beam  steering  is  the  following:  Suppose  that  the  beam  profile  is  modeled  as  a 
one-dimensional  Gaussian  (in  variable  x),  and  that  the  profile  is  displaced  a  distance 
±6.  The  average  of  such  a  profile  over  a  region  of  width  A  is  calculable  in  terms  of  the 
error  function  erf(«),  so  that  the  contrast  modulation  takes  the  analytic  form  indicated  in 
Figure  35,  for  the  special  case  that  the  averaging  region  equals  the  turbulence-limited 
resolution  cell  size.  The  contrast  modulation  due  to  beam  steering  calculated  in  this 
manner  is  plotted  in  Figures  36  and  37  for  Gaussian  beams  with  2-a  widths  of  3  m 
and  5  m,  respectively,  for  the  case  indicated  by  the  solid  line  in  Figure  32  (i.e.,  a  10- 
km  path  5  m  above  MSL)  at  a  wavelength  of  0.532  pm.  The  local-averaged  contrast 
modulation  amounts  to  ±1 4%  for  the  3-m  beam  and  ±8.4%  for  the  5-m  beam. 

Our  analysis  of  speckle  contrast  in  the  presence  of  turbulence  seems  to  be 
borne  out  by  experimental  work  done  at  the  Oregon  Graduate  Center  (Refs.  13-14) 
which  indicates  that  speckle  contrast  will  be  near  unity,  independent  of  turbulence, 
when  turbulence  is  strong  (o,2  =  o^2  »  0.1),  as  well  as  when  turbulence  is  very 
weak  (a,2  £  0.1).  In  the  intermediate  regime  of  turbulence,  Fante’s  formulation  (Ref.  - 
15)  predicts  that  the  normalized  intensity  variance  follows  the  asymptotic  trend: 

a,2  =  1  +0.99/(o,2  f* 
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Our  conclusion  is  that  for  stationary  Lambertian  targets  (such  as  hard  rocket  bodies  in 
the  visible),  the  speckle  contrast— including  turbulence  effect— will  be  >1 ,  but  for  fast- 
moving  plume  particles,  the  speckle  contrast  will  be  «1 .  These  fundamental  asser¬ 
tions  underlie  the  validity  of  the  texture  image  discrimination  concept  and  they  should 
be  capable  of  being  tested  at  ISTEF. 
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7.0  CONCLUSIONS 


From  the  analyses,  design  studies  and  experiment  plans  accomplished  so  far, 
there  appear  to  be  no  major  obstacles  to  conducting  the  experiments,  as  designed,  at 
ISTEF.  Data  from  these  experiments,  after  reduction  and  analysis,  will  allow  several 
important  goals  to  be  achieved.  First,  the  estimates  of  expected  target  returns  either 
would  be  verified  or  would  need  to  be  corrected.  Secondly,  there  will  now  be  a  data 
base  for  subsequent  high-altitude  active  tracking  experiments.  If  similar  targets  are 
used  for  balloon  experiments,  for  example,  there  will  be  a  set  of  real  numbers  to 
compare  against.  Thirdly,  it  will  be  possible  to  assess  the  feasibility  of  using  autodyne 
techniques  for  robust  active  tracking  in  the  presence  of  plumes.  By  scaling  the  data, 
one  can  predict  signal  levels,  Doppler  bandwidths,  temporal  stability  and  dynamic 
range  and  thereby  study  hardware  options  and  algorithms  to  determine  whether 
autodyne  tracking  can  be  implemented.  Fourthly,  it  is  expected  that  texture  imaging 
as  a  result  of  Doppler  bandwidth  differences  will  be  demonstrated.  These  results  will 
generate  concepts  for  new  active  discrimination  and  tracking  techniques. 

Finally,  if  Experiment  4  can  be  implemented  in  future  efforts,  it  will  serve  as  a 
test  bed  to  evaluate  interferometric  direct  detection  algorithms  against  real  images  in 
determining  position  estimate  accuracy  and  bias. 
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Figure  1.  Doppler  direct  detection  methods. 
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Figure  4.  Kinematic  implications. 
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Figure  5.  Data  quantity  (kinematic  implications). 
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Figure  7.  Experiment  1:  No  tracking,  2-m  spot. 
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*  Additions 

**  If  beam  forming  optics  will  delay  program,  omit  expanders,  use  Experiment  1  beams  with 
MBFI  in  green. 


HARDWARE 


DETECTORS 
BEAM  SPLITTERS  (1") 

FILTERS,  INTERFERENCE 
FILTERS,  GLASS 
FOCAL  PLANE  PACKAGE 
LIGHT  SHIELD,  SEPARATE  COVER 
TELESCOPE  COMPONENTS,  REVIEWED,  ON  ORDER 
WINDOW 
MIRRORS 
BARLOW  LENS 

TELESCOPE  MEASUREMENTS 

FOCAL  PLANE  ADJUSTMENTS 
BARLOW  -  FOCUS 
INCLUDE  X-Y  IN  DESIGN 


Figure  10.  Hardware  for  Experiment  1. 


Figure  11. 

DETECTOR  PACKAGES 
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Dichroic  R  =  99.7  @  0.532 

T  =  92.5  @  1 .064 


1 .  Mirror  material:  Optical  fused  silica  Homosil  or  equivalent 

2.  Spacer  material:  Zerodur,  ULE  or  equivalent 

3.  Mirrors  to  be  optically  contacted  to  spacer 

4.  Mirror  surfaces:  matched  to  better  than  Ml  20,  after  coating,  over  central  48  mm 
R  =  0.93-0.95  @  0.53  &  0.63  pm 

5.  Surface  flatness:  M20  over  central  50  mm 
R  <  0.0025  @  0.53  &  0.63  pm 

6.  Wedge  in  mirror  substrates:  5-30  arc  minutes 

Mirrors  must  be  identical  and  mounted  with  opposed  wedge  angles  such  that  the 
net  deviation  of  a  transmitted  beam  through  the  etalon  is  less  than  1 0  arc  seconds. 

7.  Spacer  thickness  to  be  measured  ±0.0025  mm 

8.  Spacer  to  have  a  wedge  of  3.3  arc  seconds  ±1 0% 

This  can  be  inspected  by  illuminating  the  etalon  with  an  expanded,  collimated  laser 
beam  @  0.633  pm.  The  fringes  seen  in  transmission  shall  be  spaced  by  20  mm 
±  2  mm. 


Figure  14.  Air-spaced  wedge  etalon. 
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Figure  16.  Transmission  fringe  sensitivity  to  incidence  angle. 


Transmission  Fringe  Sensitivity 
Nominal  Order  =  47000:  FSR  =  1 2  < 


Transmission  fringe  sensitivity  to 


o.izr 


Figure  23.  Calibration. 
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Figure  24. 

Co  PI  MfL/Tn  KO.V. 


I z"  f/tS  P.O.v/. 


Figure  26.  Calibration  board  on  gantry. 
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Range  (km) 


Figure  28. 

Signal  Photon  Counts 
Plume  and  Hard  Body 


Range  (km) 


Figure  29. 

Background  Photon  Counts 
Plum©  Fills  Receiver  FOV 


Range  (km) 


„  Figure  30. 

Lidar  Signal  to  Noise  Ratio" 

Plume  Against  Plume  Background 
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102 


Figure  31. 

EXPRESSION  FOR  SPECKLE  CONTRAST 
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Contrast  =  v/Var(/(x2))  /  <  7(x2)  > 
Triple-product  covariance  analysis  (1974) 


Figure  32. 

TURBULENCE  EFFECTS 
Spherical  Waves,  X=  0.532  jim,  Horizontal  Paths 


Range  (km) 


TURBULENCE  EFFECTS 
Spherical  Waves,  Xs  0.532  pm,  Horizontal  Paths 


Figure  33.  Predictions  for  very  low  turbulence. 


Figure  34. 

TURBULENCE  EFFECTS 
Spherical  Waves,  X=  1.064  jim,  Horizontal  Paths 
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Figure  35. 

CONTRAST  MODULATION  DUE  TO  BEAM  STEERING 

Approximate  One-DImenalonal  Analytical  Formulation 
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CONTRAST  MODULATION  DUE  TO  BEAM  STEERING 

Gaussian  a  =  1  &  m;  Ar*s  =  0.42  m 


o 

o 


an 

K 

o 


to 


to 

o 


o 

o 


uoiisinpoyi  ‘emoid  A||sueiu|  peziieuuoN 


-57- 


Position  (m) 


im  it*et 


o 


Gaussian  a  *  2.S  m;  A m  ■  0.42  m 


Table  I. 

TARGET  &  SYSTEM  PARAMETERS 


Table  II. 

BACKGROUND 


FOV  A  =  3.5x10  9  sr  <£1/2  =  250  ^trad  0  —  1.661  xlO 

_ 


APPENDIX:  RADIOMETRIC  CALCULATIONS 


Figure  A-1  defines  the  quantities  employed  in  the  radiometric  analysis.  A  far- 
field  formalism  is  employed,  i.e.,  the  laser  beam  is  represented  as  subtending  a 
transmission  solid  angle  that  is  invariant  with  range,  and  likewise  for  the  receiver 
viewing  returns  from  the  target.  The  fundamental  quantity  to  be  computed  is  the 
number  of  photoelectrons  detected  per  resolution  cell  at  the  receiver,  ;  this  is  done 
in  two  steps.  One  first  computes  the  expected  value  for  the  total  number  of  photo¬ 
electrons  in  the  entire  image  (Ntotal);  one  then  divides  by  the  number  of  resolution 
elements  (resels)  to  obtain  Nres . 

The  total  number  of  photoelectrons  is  represented  in  Figure  A-1  as  the  product 
of  four  terms:  (1)  the  number  of  photons  transmitted  by  the  laser  to  the  vicinity  of  the 
target  (including  the  outward-bound  atmospheric  loss  factor);  (2)  the  target  reflection 
cross  section;  (3)  the  fraction  of  the  scattered  radiation  which  is  collected  by  the 
receiver  entrance  pupil  (including  the  return-trip  atmospheric  loss  factor);  and  (4)  the 
receiver  photoelectron  efficiency. 

The  dependence  of  the  number  of  signal  photoelectrons  on  atmospheric 
conditions  is  summarized  in  Figure  A-2.  The  round-trip  loss  factor  is  assumed  to  be 
governed  by  the  exponential  (Beer-Bouguer)  law  shown-i.e.,  single  scattering  is 
assumed. 

The  number  of  resels  in  the  image  is  assumed  to  be  governed  primarily  by 
atmospheric  turbulence.  The  limiting  angular  resolution  0atm  imposed  by  clear-air 
turbulence  is  determined  by  Fried's  coherence  diameter  r0  through  the  approximate 
relation,  6atm  =  Ur0 . 

It  has  been  assumed  that  the  principal  source  of  noise  arises  from  the  radiance 
of  the  rocket  plume.  The  plume  is  modeled  as  a  greybody,  i.e.,  a  Planckian  radiator 
with  a  specified  emissivity  and  color  temperature.  As  indicated  in  Figure  A-3,  the  total 
number  of  background  photons  is  again  the  product  of  four  factors:  (1)  the  number  of 
photons  emitted  by  the  plume  per  unit  area;  (2)  the  observed  target  area  (product  of 
R2  and  the  solid  angle  subtended  by  the  receiver  field  of  view);  (3)  the  fraction 
intercepted  by  the  receiver  (including  a  one-way  atmospheric  loss  factor);  and  (4)  the 
receiver  efficiency.  Once  again,  to  calculate  the  number  of  photoelectrons  per  pixel, 
one  divides  the  receiver  FOV  by  the  number  of  turbulence-limited  pixels  which  it 
spans. 

Figure  A-4  indicates  the  parameter  values  employed  in  the  analysis. 
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N  total  =  Total  number  of  photoelectrons  received, 
r„  =  Fried’s  coherence  diameter  (m), 

0t  =  transmission  solid  angle  (sr)  =  tt 0 1 2/ 4 ,  and 

Nres  =  Number  of  photoelectrons  per  effective  resolution  cell. 

Figure  A-l.  Radiometric  analysis  of  signal  photocounts  per  resel . 
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Typically,  r0  should  be  several  centimeters. 
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Figure  A-4. 

RADIOMETRIC  ANALYSIS 

Parameters  Employed 

for  radiometric  analysis. 
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T|  (electrons/photon)  0.1 

r0(m)  0.05 


